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ABSTRACT 

According to the KCDM cosmological framework, galaxies underwent multiple merg- 
ers in their history. In this paper we propose to use the power spectrum of the residual 
fluctuations of the rotation curve velocity as a probe of past mergers. The proposition 
relies on the assertion that mergers are expected to induce large scale flows and in 
case of major mergers shocks are induced as well. Instabilities of the large scale flows 
and shocks could generate a large scale turbulence whose size is comparable to the 
galactic disk dimensions. We develop expressions relating underlying turbulence spec- 
tral function to the observational power spectrum of the residual of the rotation curve 
velocity. This relation can be used to test whether turbulence exists in a given galaxy. 
The method is applied to the regular spiral galaxy NGC3198 with the conclusion that 
it underwent a minor merger about 7 Gyr ago. 

Key words: galaxies: general, galaxies: kinematics and dynamics, turbulence, ISM: 
kinematics and dynamics, galaxies: mergers, galaxies: individual (NGC3198) 



1 INTRODUCTION 

The standard ACDM cosmology implies that galax- 
ies grow during their history by accretion and mergers 
(iRobertson et all 120061 : iKavirai et"afl 120091 ; IStewart et ail 
We propose to use the power spectrum of the residual 
fluctuations of the rotation curve velocity as a probe of past 
mergers. The proposition relies on the plausible assumption 
that mergers induce large scale flows, and in the case of 
major mergers shocks are induced as well. The large scale 
flows and shocks can generate a turbulence, in the ISM. The 
physical parameters of the ISM make the generation of tur- 
bulence quite easy once an instability and an energy source 
are avai lable and indeed turbulen ce in the ISM is quite ubiq- 
uitous (jElmegreen fe Scald 120041 ). The instabilities related 
to mergers are e.g. the Kelvin Helmholtz instability and the 
Richtmyer-Meshkov instability due to the shock acceleration 
|Mikaehanlll99d :l Graham fc Zhanell2000l ). 

The spatial scale of the generated turbulence is expected 
to be comparable to the galactic disk size. Such a large scale 
turbulence, is different from local turbulence generated by 
e.g. supernovae (scale of the order of 10 pc), super shells 
resulting from chains of supernovae (scale of the order of 
100 pc), or winds from early type stars (scale of the order 
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of 10 — 100 pc). Moreover, a spatial scale comparable to a 
galactic size implies for turbulent rms velocities of the order 
of 10km/ s, a turbulence life time of a few Gyr. For minor 
to medium mergers (the most common) the turbulence is 
expected to be subsonic and so its life time is even longer. 
The long life time of the turbulence implies that it can be 
detected also in nearby galaxies. Thus, the turbulence serves 
as a fossil evidence of the past merger. 

If such primordial turbulence exists, it is expected 
to be superimposed on the rotation velocity. Hence, part 
(or all) of the fluctuating residual of the observational 
rotation velocity may be identified with this primordial 
turbulence. Hints for the existence of such a turbulence are 
found in observations of r a dial ve lo cities in disks of galax 
ies (iBeauyais fc Bothunl Il999l: IWong, Blitz, fc Bosm 
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Trachternach et al.ll2008D and the observation of iBegemanl 
1 19891 ) that the residual of the rotation velocity as function 
of the azimuthal angle along a ring in the galaxy plan, 
is oscillatory. The above velocities are of the order of 
5 — Wkm/s. 

The theoretical consideration and the observational 
hints motivate one to compute the observational power spec- 
trum of the fluctuating residuals of the rotation velocity. The 
residuals depends of course on the definition of the mean ro- 
tation curve. To avoid this ambiguity we shall consider only 
galaxies with a flat rotation curve over a sizable part of the 
radial extent. The residuals are defined with respect to the 
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constant mean value of the rotation velocity. The resulting 
observational power spectrum can be used to distinguish be- 
tween random residuals (or local small size fluctuations) and 
fluctuations which are a manifestation of a galactic scale tur- 
bulence. If the latter is revealed, it means that it has been 
generated by a process with a spatial coherence scale com- 
parable to that of the galaxy. A merger or close passage by 
another galaxy is the natural candidate. 

We derive the relation between the power spectrum of 
the observational residual rotation velocity, and the under- 
lying 3D turbulence spectrum. This relation can be applied 
to determine the turbulence 3D spectral function from the 
observational power spectrum. In addition, the latter sup- 
plies an estimate for the line of sight depth of the turbulent 
region. 

The method is demonstrated for the rotation curve of 
the spiral galaxy NGC3198 and leads to the conclusion that 
in spite its ordered shape, this galaxy has experienced a 
minor merger about 7 Gyr ago. 



2 POWER SPECTRUM OF THE RESIDUAL 
ROTATION VELOCITY CURVES 

To test whether a large scale turbulence exists in a given 
galaxy one needs the relation between the power spec- 
trum of the observational residuals and the spectral func- 
tion of the underlying 3D turbulence. The observational ro- 
tation curve is usually obtained by use o f the tilted ring 
model l|Rogstad. Lockhart. fe Wright! [19741 ). We show that 
its power spectrum gives essentially the same results as a 
power spectrum of the residual Position - Velocity (PV) 
curve along a given axis on the galactic disk. The advantage 
is that the rotation velocity at a given radius, is based on 
many measurements at this radius, along the ring circumfer- 
ence and so the precision is better than that of the PV data 
which is based on a single measurement for each radius.. 



2.1 Turbulence Spectrum 

The 3D spectral function of the turbulent velocity, <E>(fc) is 
defined in terms of the 2-point correlation of the turbulent 
velocity field < v(r') -vjr + r ') > where th e angular brackets 
denote ensemble averaging (|Lesiurs!ll997h . In practice, the 
ergodic assumption is invoked and the ensemble average is 
replaced by space, area or line averages. For homogeneous 
turbulence the two point correlation is a function of the 
separation between the two points, so that 



$(fe) = 



(2tt) 3 

the inverse transform yields 



< v(r') ■ v(r + r') > e hr d A 



< v(r') ■ v(f+ r') >- 



(1) 



(2) 



In the homogeneous and isotropic case <E> depends only on 
the absolute value of the wave number and it is useful to 
introduce the turbulence energy spectrum E(k) and the tur- 
bulent velocity spectral function F(k) = 2E(k) so that 



$(fc) = $(&) = 



F(k) 

-ink' 2 



The measured radial velocity at each position in the 
plane of the sky, is an intensity- weighted average of the ve- 
locity contributed by the emitting gas along the line of sight. 
A simplifying assumption of homogeneity translates this into 
a line of sight average of the velocity. 

In the following two subsections we assume that a 3D 
large scale turbulence exists and derive the relations between 
the observational power spectrum of the PV and rotation 
curve data and the underlying 3D turbulence. The simplify- 
ing assumptions of homogeneity and isotropy are used. 



2.2 The Power Spectrum of Position Velocity 
Data 

Let us consider first the power spectrum of a Position- 
Velocity (PV) data, namely the rotation velocity as function 
of position along a given axis of the galaxy disk. This power 
spectrum is defined as Fourier transform of the two point 
correlation of the observed residual velocity. Denoting the 
axis as y = and the line of sight as z, and the residual 
observed rotation velocity as u, the power spectrum is 



Ppv(q) = — / e rqx < u(x + x',0,z)u(x',0,z') > dzdz' dx(A) 
2tt J 

Using the isotropy assumption and Eq. [2] one finds 
<u(x+x',0,z)u(x',0,z) >=\J *(fc)e- ifea;+fei(z - z ' )) d 3 fc(5) 
which when inserted in Eq. [3] leads to 

Ppv(q) = — / <&{q,k y ,k z )e kz{ - z ~ z ' ' dk y dk z dzdz (6) 



and Integrations over z, z yield 



i(fc z g/2) V 



Ppv(q,D) = lj $(q,k y ,k z ) i^^LZLj dk y dk z (7) 

with D, the depth along the line of sight of the turbulent 
region. In the case of a disk galaxy D cos i = 2H with H 
denoting the scale height of the disk and i the inclination 
angle. 

For an homogeneous and isotropic turbulence, the 3D 
turbulence spectral function is a power law of the absolute 
value of the wavenumber: $((/, k y , k z ) = A(q 2 + k y + k z )~ m . 

In the incompressible subsonic case the s pectrum 
is th e Kolmogorov spectrum with m = 11/6 l|Lesiurs! 
Il997[). For compressible supersonic turbulence m = 
2 dPassot. Pouauet. fc Woodward] 1 19881 ; iGirimaii fc Zhoul 
Il995l) . The steeper slope is due to the fact that a fraction 
of the turbulent kinetic energy density at a given wavenum- 
ber is converted to compression work decreasing the rate of 
energy transfer to the larger wavenumbers. 

Integration over k y from — oo to oo yields 



P pv (q,D) = B(m) / (q z +k 



,2s(-m+l/2) 



an(fc z -D/2) V 
(k z D/2) ) 



dk z 



where B(m) is a constant depending on m and D. One can 
further express the power spectrum in the form 



k = Ifcl 



(3) 



P pv {q,D) = B{m) (qD/2) 2 + V 



n 2^(-m + l/2) SHI 7] 

r/ 2 



d V (8) 
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The power spectrum is a power law with index 2 — 2m for 
qD/2 « 1 and 1 - 2m for qD/2 » 1. This effect of the 
change of the power law index was found o bservationally by 
lElmegreen. Kim, fc Staveley^S mith (2001) for the Hi power 



spectrum of the LMC, by iDutta et all (|2009l ) for the Hi 



power spectrum of NGC1058 and recently by iBlock et al.l 
(2010 ) for the infrar e d pow er spectrum of the LMC, and by 
IContini fc Goldman! l|2010l ) for both the velocity and mid 
infrared power spectra of the shocked nebulae near the tur- 
bulent Galactic Center. The results of the last two references 
are based on the Spitzer data. 



2.3 The Power Spectrum of the Rotation Curve 

The rotation curve is obtained by use of the tilted ring model 
(jRogstad. Lockhart. fc Wrightll 19741 ) in which the disk is ap- 
proximated as a superposition of concentric rings or annuli. 
The observed velocity field is then used to derive for each 
ring a rotation velocity, an inclination angle and a position 
angle. The rotation curve velocity as function of radius ob- 
tained in this way is more accurate than the PV curve, since 
the value at each radius is derived from multiple observa- 
tional values. 

This motivates a computeation of the power spectrum 
of the residual fluctuations of the rotation curve and its use 
to test for the existence of an underlying 3D turbulence. 
For random fluctuations the power spectrum is expected to 
be independent of the wavenumber. If however, the power 
spectrum varies in an ordered manner with the wavenumber 
it is possible that it reflects an underlying turbulence. 

As in the case of the PV power spectrum there is the 
need to obtain the relation between the observationally com- 
puted power spectrum and the underlying turbulence spec- 
tral function. In the case of the rotation curve this task 
seems a prioiri more difficult than the PV case since the 
rotation velocity is fitted to the observed line of sight velci- 
ties along the ring circumference and not measured durectly 
along a given line in the plane of the disk. For simplicity the 
fit along the ring is represented as a simple average over its 
circumference, so that 



P{9) 



(2n) 



-} qr < u(r+r',8,z)u(r',8',z') > dOdO 1 dzdz (9) 



As in the PV case, here too use of Eq. [2] and assuming 
isotropy of the turbulence imply 



3 (2tt) 3 



$(fc k z )e~ ik< - T+r '^ cos ( a_e ) e ikr ' cos ( a - e ')^Q^ 



iqr —ikr(z- 



-*') 



drdakdkdk z dzdz d9d6 



where k = yfcl + fcf is the wavenuber in the plane of the 
sky and a denotes its azimuthal angle. In the integration 
along the ring circumference a spans an interval of 2n and 
so do a — 8 and a — 6' . One gets, 
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Figure 1. P(q) in dimensionless units, as function of qD/2 for 
m = 11/6 



$(k,k z )e iqr J (kr)k 



■m(fc.£>/2) \ 3 
(k z D/2) ) dkdTdkz 



with Jo denoting the zero order Bessel function of the first 
kind. 

For an underlying 3D turbulence characterized by a 
power law spectral function &(k,k z ) — A(k 2 + k 2 )~ m , the 
integrations on r and k yield 

, , , 2 ,(- m+1/2 ) ( sin(k D/2) \ 2 dkz{l2) 
(k z D/2) J 



P(q) = B^m) 



where Bi(m) is a constant depending on the index m and 
D is the depth along the line of sight. The analytic integra- 
tions over r and k le ading to Eq. ()12[) were done using the 
[Mathematical (j2010h software. One can further express the 
power spectrum in the form 



P(q) = B 1 (r 



(qD/2) 2 + V 2 ) ( ~ m+1/2) 



sin 2 r\ 
V 2 



dr/ (13) 



Note that the observational power spectrum of the residual 
velocity curve has the same relation with respect to the un- 
derlying 3D turbulence as does the observational PV power 
spectrum. This result should not come as a surprise since 
the rotation curve can be regarded as an averaging of PV 
data along many axes. As a result, the power spectrum of its 
residual fluctuations is an average of the PV power spectra 
and hence share the same relation with the spectrum of the 
underlying 3D turbulence. As noted above, it is expected to 
be more accurate than the PV power spectrum. 

Fig. Q] displays the observational power spectrum as 
function of qD/2 for an underlying 3D Kolmogorov turbu- 
lence characterized by m = 11/6. For spatial scales much 
larger than the depth along the line of sight, the index of the 
power spectrum is —5/3 and for spatial scales much smaller 
than the depth it approaches —8/3. The transition between 
the two regimes is at kD/2 ~ 1.4 and there the logarithmic 
slope is —13/6. 



P(q) 



3 (2tt) 3 



§(k,k z y qr e lk * (z - z ' ) 



(11) 



lkrcos f : ' drd/3kdkdk z dzdz' ' d9d9' 



2.4 The Observational Power Spectrum of the 
Rotation Curve of NGC9183 

The proposed method is demonstrated by applying it to 
NGC9183. The latter appears to be a rather ordered 
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Figure 2. The residual velocity in units of km/s as function of 
radial distance in units of arcminutes. The error bars are taken 
from Table 2. of Begeman (1989). 

galaxy, although some deviations from regularity were noted 
l|Bosma|[T98ll ; Ivan Albada et al. l ll985l : lBegemarilll989r i. The 
Hi column densities do not exceed 5 x 10 20 cm -2 so the 
optical d epth is ra t her s mall. The data presented in Ta- 
ble 2. of iBeeemanl (|l989h show a very flat rotation curve 
for radial distances in the range 1.5 — 11' corresponding to 
4.1 kpc — 29.9 kpc. This is the range of radial distances to 
which the proposed test is applied. 

The residuals (with respect to the mean value) are plot- 
ted in Fig. 2 as function of radial distance in units of ar- 
cminutes. Up to 3', the radial distances are spaced by 0.25' 
and later by 0.5'. Since the computation of power spectra re- 
quires even spacings, we omitted the data points at positions 
1.75', 2.25', 2.75' thus there remain overall 20 evenly spaced 
data points. Note that the uncertainties are (2 — 3)fcm/s; 
substantially smaller than the 15.5 km/s uncertainty in a 
single measurement l|Begemanl Il989l ). This arises from the 
use of many measurements along the ring and demonstrates 
the advantage of using a rotation curve rather than a PV 
data. 

The residual velocity shown in Fig. [2] does seem to ex- 
hibit fluctuations, as function of radial distance. But one has 
yet to test whether these fluctuations do indeed represent a 
large sca le turbulence. 

The lMathematica7l (|2010h software was used to com- 
pute the power spectrum of the residual rotation curve ve- 
locity as function of the wavenumberwas q — 2n/l, with 
I denoting the corresponding spatial scale. Because of the 
cyclic nature of the discrete Fourier transform, the power 
spectrum is obtained for relative wavenumbers in the range 
1 — 11 with q = 1 corresponding lo = 25.8fepc, the largest 
spatial scale of the turbulence and q — 11 corresponds to 
2.35 kpc, which is about twice the radial spacing between 
the observations used in the present analysis. 

To obtain the error bars for the observational power 
spectrum we performed 10 4 simulations of "observational" 
velocity sets. At each position a velocity was randomly cho- 
sen from a normal distribution with a mean equal to the 
observational value at this position and a standard devia- 
tion equal to the observational uncertainty at this position. 

For each such set, the power spectrum was computed 
and subsequently the standard deviations, at each wavenum- 
ber, of the logarithm of the simulated power spectra. These 
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Figure 3. Dots:Thc logarithm of the power spectrum of the 
residuals of the observational rotation curve in units of (km/s) 2 , 
as function of the logarithm of the relative wavenumber q. lower 
curve: the best fit <I>(g): has depth D = 7.3lkpc corresponding 
to an Hi scale height of 1130 pc. upper curve: the best fit &(q) 
without the q = 7 point has depth D = 3.78kpc corresponding 
to an Hi scale height of 584 pc. 



standard deviations were taken as the uncertainties in the 
logarithm of the observational power spectrum of the resid- 
ual rotation velocity. 

The standard deviation of the observational residual 
rotation velocity of Fig. [2] is 3.6km/s. Hence, if it indeed 
represents a turbulence the latter must be subsonic(in the 
isotropic case Vturb = \/33.6 = 6.2 km/s). The theoretical 
fit is thus chosen to be a turbulence spectral function with 
m = 11/6. The depth D, along the line of sight and the 
overall normalization of the function given by Eq. (|13|) to 
minimize the reduced x 2 value. 

The observational power spectrum of the residual ro- 
tation curve velocity is shown in Fig. [3] together with two 
such fits. The lower one is the best fit with x 2 = 0.43. It 
has a depth along the line of sight D = 7.31 kpc translated 
into a scale height for neutral hydrogen Hhi = 11 30 pc. The 
law va lue of \ 2 suggests that the uncertainties in iBegemanl 
(1989) are probably overestimated. 

This analysis implies that the residuals are consistent 
with being a manifestation of a subsonic turbulence with a 
largest scale of 25.8 pc, a turbulent velocity of 3.6km/s, and 
a lifetime scale of the order of IGyr. 

About 75% of the x' 2 value is contributed by the point 
at relative wavenumber q = 7 corresponding to a scale of 
about 3.7 kpc. The upper curve in Fig. [3] is the best fit 
to the power spectrum with this point excluded. It has a 
X 2 = 0.11, and a depth D — 3.78 kpc corresponding to a 
scale height Hhj = 584 pc. This latter is more consistent 
with observational values for disk galaxies (Bagctakos et al] 
120101 ). For each of the two fits there is a correlation between 
the normalization of the fit and the value of D, resulting 
in a flat dependence of x 2 on D. Thus, in each case the 
calculated D and the corresponding scale height vary in a 
range of about ±20%. 
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3 DISCUSSION 

This paper puts forward the idea that analyzing the residual 
of the rotation curve velocities could test for the existence 
of a galactic-scale turbulence. If such a turbulence exits, the 
large spatial scale implies a lifetime of the order of a few 
Gyr. These two features point to a primordial merger as 
the natural generating mechanism. The test relies on the 
fact that random residuals of the velocity curve would have 
produce a flat power spectrum in sharp contrast with known 
turbulence spectra that exhibit a power law dependence on 
wavenumber, for wide ranges of the latter. 

In this paper the method is proposed for galaxies ex- 
hibiting an observational rotation curve which is flat over siz- 
able range of radial distances. This limitation avoids ambi- 
guity as to the definition of what constitutes the residual ve- 
locity curve. In any case, if the derived observational power 
spectrum can be identified with turbulence spectrum, the 
issue of the definition of the residual is a posteriori solved. 

We derived a semi-analytic relation between the obser- 
vational power spectrum and the spectral function of the 
underlying 3D turbulence. We proved that the relation is of 
the same form for the observational spectra calculated from 
a position velocity data and from the rotation curve. The 
latter is of course advantageous since at each radial distance 
it represents a fit along a ring and not a single value. The 
method is capable of determining the depth along the line 
of sight of turbulent gas and hence the disk scale height. 

The method has been demonstrated for the rotation 
curve of NGC3198. The computed power spectrum fits very 
well a Kolmogorov spectrum with a logarithmic slope of 
—5/3 for the largest scales and changing to —8/3 for the 
small scales. The turbulent rms velocity is about 6 km/s; 
the lifetime is of the order of 7 Gyr and the implied scale 
height about 580 — 1100) pc. For this galaxy, the obser- 
vationa power spectrum spans a decade in wavenumbers. 
Clearaly, an higher spatial resolution would have enabled a 
wider wavenumbre range, inmprove the numerical fit of the 
the depth and thescale height, and lend more credibility to 
the cocnclusions. But even so, the case for an underlying 3D 
turbulence in NGC3198 seems quite good. 

The method can be used as an observational tool for 
finding how common were mergers. It easy to identify merges 
in cases of galaxies with a disturbed shape. These mergers 
are usually major mergers. The majority of the cosmological 
mergers are expected to be minor mergers that don't have 
a sharp morphological imprint. The present method enables 
detecting turbulence also for minor mergers and for galaxies 
that appear regular, as in the case of NGC3198. The long 
lifetime of the turbulence enables using it as a fossil evidence 
for primordial mergers. 

It would be interesting to look for turbulence in ad- 
ditional galaxies and specifically, to test for a correlation 
between the disk thickness and the turbulent rms velocity. 
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